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ABSTRACT: Computing earthwork volumes is a necessary activity for nearly all construction projects and is often 
accomplished as a part of route smreying, especially for roads. Dredging and highways. If the area of the cross section at 
each station is known, then multiply that average end area by the known horizontal distance between the stations to 
determine the volume of cut. We present the results of precise static and kinematic point-positioning solutions obtained in 
post-processing, with a recursive Kalman filter estimator. The observations are dual frequency CPS carrier phase and 
pseudo-range, treated as two distinct data types, each with its own measurement equations. The receiver data is double 
differenced between satellites, to eliminate the receiver clock, and processed using precise satellite ephemerides and clock 
corrections. While estimating receiver coordinates, ionosphere-free carrier phase Lc biases, and tropospheric refraction 
model errors were treated accordingly. The motivation has been to show the possibility of the application of GPS in 
Topographical Survey and it associated accuracy. The method is fast, accurate and data was acquired using Trimble 
4000SSI GPS Receiver and processed using Trimble Geomatics office.. The results indicate a precision at the level of a few 
centimeters, for static solutions, and below one decimeter, in kinematic mode. 
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I. INTRODUCTION 

Society recognises the need to move towards sustainable development. Construction has an important role to play 
within the sustainable development agenda not only because of the contribution to the national economy, but also because 
the built environment has a major impact on the quality of all of our lives, our comfort and security, our health and wellbeing 
and our productivity. 

Road construction, maintaining and upgrading our existing Road has a potential environmental impact. The 
challenge today is to deliver road that can stand test of time and enhance the quality of life, while at the same time reducing 
the suffering encountered while embarking on a journey. 

The purpose of the paper is to discuss the Geomatics modern approach used in the construction industries. 

Earthwork computations are based on the design cross sections generated in the design process and cut and fill end 
areas computed and stored by that process. Earthwork computations are based on the design cross sections generated in the 
design process and cut and fill end areas computed and stored by that process. Additional required and optional data is 
needed for the earthwork computations. When subsurface data is provided, end areas for each material in cut are computed 
and is used as a basis for volumes. A compaction factor may be specified for each material, and volumes and mass ordinates 
will be computed accordingly. The earthwork parameters, such as compaction factors, added quantities, and forced balance 
stations, are used in computing the earthwork quantities for the designed roadways. These parameters are entered along with 
the other design data. These parameters are necessary in the compilation of earthwork volumes, mass ordinates and haul 
computations. Although these data are not strictly design data, they are entered and stored along with design data because the 
designer may then want to modify volume computation parameters and make new alternate volume computations. 

A compaction factor is a characteristic assigned to a specific earth material which indicates how, and to what 
extent, its volume will change when that material, after having been cut from the original ground, is compacted as fill 
material. This characteristic is expressed as the ratio of a given volume of this material, in its original ground state, divided 
by the volume that same material would occupy when compacted as fill. Compaction factors produce an adjusted fill or cut 
volume to be used in computing earthwork volumes and mass ordinates. The compaction factors are applied when 
computing volumes. Volumes are computed separately for each roadway, and the appropriate compaction factor is selected 
and applied. Then these individual roadway volumes are summed to compute the design cross section volume. 

II. GLOBAL POSITIONING SYSTEM (GPS) 

Navigation and positioning is the art of determining position, speed and orientation of an object. The GPS program 
me was born in 1972 out of the quest for the US Navy and Air force to develop all weather Global Radio Navigation System. 
The GPS was designed to be a passive survivable continuous system which will provide any suitably equipped user with 
three dimensional position (x, y, z) velocity and precise time information. 

The principle involves the measurement of distance (or range) to at least three Satellites whose X, Y and Z position 
is known, in order to define the user's X p , T p and Z p position. In its simplest form, the satellite transmits a signal on which the 
time of its departure (fD ) from the satellite is modulated. The receiver in turn notes the time of arrival (fA) of this time mark. 
Then the time which it took the signal to go from satellite to receiver is (fA - fD ) = At (called the delay time). The measured 
range R is obtained from (Andrew et.al 2004). 
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Rl = (tA-tD) C = AtC (1) 

where c = the velocity of light. 

Whilst the above describes the basic principle of range measurement, to achieve it one would require the receiver to 
have a clock as accurate as the satellite's and perfectly synchronized with it. As this would render the receiver impossibly 
expensive, a correlation procedure, using the pseudorandom binary codes (P or S), usually 'S', is adopted. The signal from 
the satellite arrives at the receiver and triggers the receiver to commence generating the S-code. The receiver-generated code 
is cross-correlated with the satellite code. The ground receiver is then able to determine the time delay (At) since it generated 
the same portion of the code received from the satellite (Abdel -Salem 2005). 

However, whilst this eliminates the problem of an expensive receiver clock, it does not eliminate the problem of 
exact synchronization of the two clocks. Thus, the time difference between the two clocks, termed clock bias, results in an 
incorrect assessment of At. The distances computed are therefore called 'pseudo-ranges'. The effect of clock bias, however, 
can be eliminated by the use of four satellites rather than three. A line in space is defined by its difference in coordinates in 
an X, Y and Z system: 

R = (AX 2 + AY 2 + AZ 2 )" 2 (2) 
If the error in R, due to clock bias, is (_R) and is constant throughout, then: 
Ri+_R= [(Xi - X p ) 2 + (Yi- Y p ) 2 + (Z x - Z p f] 1/2 
R 2 + _R= [(X 2 - X p ) 2 + (Y 2 - Y p ) 2 + (Z 2 - Z p ) 2 ] 1/2 

R 3 + _R= [(X 3 - X p ) 2 + (Y 3 - Fp) 2 + (Z 3 - Z p ) 2 ] 1/2 (3) 
R 4 + _R = [(X 4 - X p ) 2 + (Y 4 - Y p ) 2 + (Z 4 - Z p ) 2 ] 1/2 

Where X n , Y n , Z„ = the coordinates of satellites 1, 2, 3 and 4 (n = 1 to 4) 
X p , Y p , Z p = the coordinates required for point P 
R„ = the measured ranges to the satellites 

Solving the four equations for the four unknowns X p , Y p , Z p and _R eliminates the error due to clock bias. Whilst the 
use of pseudo-range is sufficient for navigational purposes and constitutes the fundamental approach for which the system 
was designed, a much more accurate measurement of range is required for positioning in Geodesy. This is done by 
measuring phase difference by means of the carrier wave in a manner analogous to EDM measurement (Schofield 1993). 

III. KINEMATIC POSITIONING ALGORITHM 

Kinematic positioning and alignment relies on the relationship of the earner phase observations to the range 
observations described in the following equation (Raquet 1998). 

PR X =R + bu { + bsv PRl + T + I l + nPm 
f y R 2 = R -+- bu 2 -+- bsv PR2 -hT + I 1 ^f- +„ PR2 

- CPH l =R + bu CPH1 + bsv CPHX +T-I 1 + nCPm -N^ (4) 

— CPU , R -+- bu cm2 -+- £>xv cm2 ^% + nCTW2 — 7V 2 >L, 

Where 

PR = pseudo-range on LI or L2 frequencies (meters), CPH = carrier phase on LI or L2 frequencies, (meters), RT = 
true range (meters), bu = range equivalent receiver clock offset (meters), bsv = range equivalent satellite clock offset 
(meters) 

T = tropospheric delay (meters), I = ionospheric delay (meters), n = measurement noise (meters), N = CPH integer (cycles), 
X = carrier wavelength (meters) 

The ionospheric delay is different on the LI and L2 observations as it is inversely proportional to the frequency 
squared and so can be removed from the PR by differencing. The DGPS corrections will remove any errors in the navigation 
solution caused by satellite position and clock offsets. The accuracy of the PR derived DGPS corrected position solution is a 
function of the pseudo-range noise, which includes receiver noise and multipath errors. The GPS/inertial navigation solution 
will filter the short-term noise effects, but it cannot correct for correlated noise errors from multipath. 

The major problem with static GPS is the time required for an appreciable change in the satellite/receiver geometry 
so that the initial integer ambiguities can be resolved. However, if the integer ambiguities could be resolved (and constrained 
in a least squares solution) prior to the survey, then a single epoch of data would be sufficient to obtain relative positioning 
to sub-centimeter accuracy. This concept is the basis of kinematic surveying. It can be seen from this that, if the integer 
ambiguities are resolved initially and quickly, it will be necessary to keep lock on these satellites while moving the antenna 
(Schofield 1993). 

IV. RESOLVING THE INTEGER AMBIGUITIES 

The process of resolving the integer ambiguities is called 'initialization' and may be done by setting-up both 
receivers at each end of a baseline whose coordinates are accurately known. In subsequent data processing, the coordinates 
are held fixed and the integers determined using only a single epoch of data. These values are now held fixed throughout the 
duration of the survey and coordinates estimated every epoch, provided there are no cycle slips. The initial baseline may 
comprise points of known coordinates fixed from previous surveys, by static GPS just prior to the survey, or by 
transformation of points in a local coordinate system to WGS 84 (Abdel - Salam 2005). 
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An alternative approach is called the 'antenna swap' method. An antenna is placed at each end of a short base (5-10 
m) and observations taken over a short period of time. The antennae are interchanged, lock maintained, and observations 
continued. This results is a massive change in the relative receiver/satellite geometry and, consequently, rapid determination 
of the integers. The antennae are returned to their original position prior to the surveys. It should be realized that the whole 
survey will be invalidated if a cycle slip occurs. Thus, reconnaissance of the area is still of vital importance, otherwise re- 
initialization will be necessary. A further help in this matter is to observe to many more satellites than the minimum four 
required (Andrew et.al 2004). 

IV. 1. KALMAN FILTER 

The Kalman filter is a digital algorithm that provides current estimates of a system variable such as position 
coordinates (Schoffield, 1993). The filter uses statistical models to properly weigh each new measurement relative to past 
information. It also determines up - to - date uncertainties of the estimates for real - time quality assessments (Kaplan 
1996). 

The Kalman filter is a multiple input multiple - output digital filter that can optimally estimate in real time the 
states of a system based on its outputs. These states are all the variables needed to completely describe the system behaviour 
as a function of time such as position and velocity. The Kalman filter filters noisy measurements to estimate the desired 
signals. The estimates are statistically optimal in the sense that they minimize the mean square estimation error. The Kalman 
filter contain a dynamic model of the GPS receiver platform motion and outputs a set of user receiver position, velocity and 
Time (PVT) state estimates as well as associated error variances [6] 

The dynamic model can be derived by a Taylor's series expansion about the true position of the receiver. We let u 
(t) represent the true position of the receiver at a time t. Then at a time t; shortly after a time tu the receiver will be at a 
position (Kaplan 1996). 

tt, \ tt, x du(t)\ . . 1 du 2 (tj , , 2 1 dw'im , , 3 

U{t) = U{t 0 ) + —^A (t-t Q ) + — M (t-t 0 ) + -4f (t-t 0 ) (5) 

^T|r=, 0 2! dt 2 |, =to 0 3! dt 3 \, =lo ° 

Where 

—^A = Velocity (5a) 
dt \,= tv 

1 du(t) 

7T~ = acceleratt on (5b) 

2! dt 2 ,=,„ 



1 du\t) . , 

The third term in the expansion (Jerk) is regarded as negligible. The non negligible terms depend on the system 
being modeled and may be estimated for position measurement. Filters designed for PVT determination typically estimate 



either user states position (x u , y u , z u ) velocity X u> V u Z u receiver clock offset (tu) and receiver clock drift 



IV.2. CODE MEASUREMENT 

The GPS code measurement p is a linear measurement known as pseudo range in meters. Pseudo Range is a 
measure of the distance between the satellite and receiver at the time epoch of transmission and reception of the signals. 

The transmit time of the signals is measured by comparing (correlating) identical Pseudo Random Noise (PRN) 
codes generated by the satellite and the receiver (Ehiorobo 2009). 

p = R + C (dt - dT) + d lon + d t rop + d mpath . (6) 

where 

R = True range between the S V (at transmit time) and receiver (at r receiver time) 



C = Speed of light (m/s) 

dt = Receiver clock error (sec) 

dT = Satellite clock error (sec) 

d ion = Ionospheric delay parameter 

dt rop = Measurement delay due to troposphere (m) 



dmpath = Measurement delay due to multipath 

Most DGPS Processing involves generating a nominal "computed" range between the receiver and the satellite 
which is calculated from the known coordinate of the receiver locations and the satellite [6]. 

The true range vector R is 
R = P sv - P rec " (7) 

Where; P sv - is the computed position of the satellite at transmit time (From ephemeris) 
P rec is the nominal position of the receiver at receiver time. 
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IV.3. CARRIER MEASUREMENT 

The constant motion of the GPS satellite constellation in orbit demands that the receiver be capable of taking care of 
the changing Doppler frequency shift on Li. 

In the case of Dual frequency receivers both Lj and L 2 are tracked. L! wavelength is 19cm and L 2 wavelength is 
24cm. The shift in frequency arises due to the relative motion between the satellite and the receivers. Integration of the 
Doppler frequency offset results in an extremely accurate measurement of the advance in signal carrier phase between time 
epochs (Gao et.al 2001). 

Using Double Differencing Processing DD techniques on C/A or P(Y) - code earner phase observables removes 
most of the error sources. 

To achieve centimeter level accuracy, it is still necessary to carry out further refinement of the propagation path length 
measurement using carrier - cycle integer ambiguity resolution. 

Once the receiver locks on to a particular satellite, it not only make C/A and or p(Y) code pseudo range 
measurements on Li and L 2 (if L 2 is present) but it also keeps a running count based upon the Doppler frequency shift 
present on Li and L 2 carrier frequencies. 

Each epoch for running this cycle count is available for the receiver. It should be noted that the advance in carrier 
phase during an epoch is determined by integrating the carrier Doppler Frequency offset (fo) over the interval of the epoch. 
Frequency fo is the time rate of change of the carrier phase, hence integration over an epoch yields, the carrier phase advance 
during an epoch. At the conclusion of each epoch, a fractional phase measurement is made by the receiver. This 
measurement is derived from the carrier phase tracking loop of the receiver. Mathematically the relationship is as follows 
(Abdel - Salam): 

= = j fo,,(r)dr + 0 IAO 

't (8) 
Where fa lo = O 

Where fa 20 = O 

Where 

(j) is the accumulated phase at the epoch shown 
n and n-1 are the current and immediate past epochs 
f D - is the Doppler frequency as a function of time 
§i is the fractional phase measured at the epoch shown 



IV.4. IONOSPHERIC DELAY 

The ionosphere is the part of the upper atmosphere where free electrons occur in sufficient density to have an 
appreciable influence on the propagation of radio frequency electromagnetic waves. It is a weak ionized plasma extending 
from approximately 50Km to lOOOKm above the surface of the earth and can affect radio wave propagation in various ways 
(Kaplan 1996). 

Ionosphere is a dispersive medium. The refractive index is a function of the operating frequency. The refractive 
index n of the ionosphere can be expressed as (Schofield 1993). 

X 



2 =1- 



2(1- X -iZ ) 



[4(1- x-iZf J 



Cos 9 

' f 



Z = — ;w = ZTLf 

and 

f - is frequency of incoming signal in Hz 

N e - is the electron density in electron / m 3 

e - is the electron change = -1.602 X 10" 19 Coulomb 

e n - is the permeativity of free space = -8.854 X 10" 12 Farad / m 

m - is the rest mass of electron = 9. 107 X 10~ 31 Kg 

6 - is the angle of the ray with respect to the earth's magnetic field. 

V - is the electron - neutron collision frequency and 

f H - is the electron gyro frequency typically 1.5MHz 

IV.5. TROPOSPHERIC DELAY 

Tropospheric delay occurs as a result of the GPS signals passage through the Earth's lower atmosphere. This 
region extends to an altitude approximately 10km above the earth's surface. The stratosphere is also considered part of the 



www.iimer.com 



2290 I Page 



International Journal of Modern Engineering Research (IJMER) 
www.iimer.com Vol. 3, Issue. 4, Jul. - Aug. 2013 pp-2287-2295 ISSN: 2249-6645 

lower atmosphere. It extends upward from the tropopause to an altitude of approximately 50 kilometers in these two sections 
of the atmosphere, electromagnetic refraction slows the GPS radio signals causing a delay in its arrival at the receiver [6]. 

IV.6. RECEIVER NOISE 

Receiver noise can be considered as white noise as it is uncorrelated over time. In the case of GPS, it is the error in 
phase and code measurement due to imperfect tracking of the GPS signals by the phase and delay lock loops. The tracking 
error is random and is assumed to have a Gaussian distribution. The magnitude of the noise is a function of the measurement 
type and the signal strength. 

The C/A - code measurement can be tracked with a typical accuracy of 1.5 meters but it may vary between 0.2 and 
3 meter depending upon the strength of the receiver signal. Typical receiver noise for P(Y) code measurement is 10 - 30cm. 
Code and carrier phase measurement noise can be smoothed by an estimate of a "zero baseline" test (Kaplan 1996). 
The GPS signal is split into two and fed to two separate but same types of receivers. In static application, the predominant 
error component in earner phase noise is jitter of the phase lock loop caused by thermal noise. This can be expressed 
according to; (Schofield 1993). 




1 

2TC 



(10) 



where 

B n - carrier loop noise band width in Hz 

c c/No/10 

— - is the carrier to noise power expressed as 10 with C/No in 

n i0 

dB- Hz 

T - is the predetections integration time in seconds 

X - is earner wavelength 

IV.7. SATELLITE ORBIT ERRORS 

Satellite orbit errors result from the uncertainty in orbital information. To compute the satellite position, we either 
use broadcast ephemeris or precise ephemeris. Although precise ephemeris information is much more accurate than those of 
broadcast information. The satellite positions are computed from a set of Keplerian orbit and perturbation parameters with 
clock parameters which are predicted states for the satellite orbit and are updated every two hours (Raquet 1998). 

IV.8. MULTI PATH 

Multi path refers to the existence of signals reflected from objects in the vicinity of a receiver's antenna that corrupt 
the direct line-of-sight signals from the GPS satellites, thus degrading the accuracy of both code -based and carrier phase- 
based measurements. The object may be tall buildings, television and Telephone masks, etc 

IV.9. GEOMETRIC DILUTION OF PRECISION 

In GPS observations, the position of the satellite affects the three dimensional angles of intersection between the 
satellite and the users position. When the satellites are close together or in a straight line, a low accuracy fix is obtained. 
When they are wide apart almost forming a square, a high accuracy is obtained (Abdel-Salam, 2005) 

Thus, the dilution of precision is used to measure user position accuracy. The satellite configuration geometry with 
respect to the ground station is called the GDOP the geometric dilution of precision defined as (Ehioro 2008); 



gdop = — Ja x 2 + (J ' + a ' + (in 

(7 " 

Where; o is the mean square error of the pseudo range, which has a zero mean 
o x o y o z - measured r ms errors of the users position in the xyz- direction and 
o b is measured rms user clock error expressed in distance. 
Other DOP parameters include 

VDOP - vertical Dilution of Precision - One dimensional (1-D) 

HDOP - Horizontal Dilution of Precision - Two dimensional (2-D) 
PDOP - Position Dilution of Precision - Three dimensional (3-D) 
The position dilution of precision is defined as 

PDOP= I (12) 

CT V 

The Horizontal Dilution of precision is defined as 

HDOP= +ct; - (13) 
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The Vertical Dilution of precision 

VDOP = ~- ^ 

CT 

The time dilution of precision 

TDOP = °~ h ^ 15 ^ 

CT 

The smallest dilution of precision DOP value means the best satellite geometry for calculating user position. 
The observation is at the base of the tetrahedron. Under this condition, the DOP values are (Raquet 1998) 
GDOP= A /3 _ « 1.73 
PD0P = 2. ' 
HDOP = VDOP 




X V3 

And finally TDOP = \ 



V. CYCLE SLIPS 

A cycle slip is a break in continuous satellite tracking. The wave fronts that the receiver counts during continuous 
tracking are also referred to as cycles. Therefore, when the receiver loses count because of an interruption in satellite 
tracking, the cycle or wave counting slips. 



V.l RESULTS OF ADJUSTMENTS OF OBSERVATIONS 

The qualities to be obtained after field measurements and adjustment in a GPS network include baseline vector 
components and their covariance. 

The final outcomes of baseline adjustment are the adjusted baseline vectors, estimated points coordinates and their 
covariance matrices. We present below results of adjustment using Trimble Geomatics office. 

Table 1 - Baseline Processing Summary 



Baseline Summary B120 (3842 to AE000091) 

Processed: 

Solution type: 

Solution acceptability: 

Ephemeris used: 

Met Data: 

Baseline slope distance: 
Elevation mask: 
Variance ratio: 
Reference variance: 
RMS: 

Horizontal Precision 1-sigma (scaled): 
Vertical Precision 1-sigma (scaled): 



Processing Summary 



Thur-da\. Jul 21. 201 1 07:44:35 P.\l 

LI fix 

Multiple 

Broadcast 

Standard 

2538.933m 

1 3 decree-. 

0.0 

I 1.688 
0.008m 
0.030m 
0.068m 



Table 2 - Baseline Components (Mark to Mark) 



From: 


|3842 


1 1 


1 1 






(Grid 


1 


|Local 


|WGS 84 






Northing 


|367124.091m 


|Latitude |7°19'01.50193"N 


Latitude 7°19'02.82741"N 




Easting 


|436929.511m 


|Longitude |6°22'02.87020"E 


|Longitude |6°22'00.44966"E 




Elevation 


|365.672m 


|Height |365.672m 


Height 366.558m 






To: 


|AE000091 










(Grid 




|Local 


|WGS 84 






Northing 


|368263.699m 


|Latitude 7°19'38.89654"N |Latitude 


|7°1940.21925"N 


Easting 


|434660.811m 


Longitude |6°2049.0781 


8"E |Longitude 


|6°2046.65650"E 


Elevation 


|308.438m 


|Height |308.438m 


(Height 


309.357m 




Baseline: 


U Northing 


1 1139.607m 


|NS Fwd Azimuth |296°54'33" 


|n x 


|48.700m 




U Easting 


|-2268.701m 


|EU. Distance |2538.153m 


|n y 


|-2272.047m 


U Elevation 


|-57.235m 


!□ Height |-57.201m 


|n z 


1 1132.082m 



Table 3 - Standard Errors 
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Baseline Errors: 


□ m Northing 


0.005m 


□ NSFwd 
Azimuth 


|o.517 seconds 


□ 


: 


X 


0.020m 


|n CD Easting 


|0.009m 


|D EU.Distance 


0.008m 




: 


Y 


|0.008m 


!□ CD Elevation 


|0.020m 


|D CD Height 


|0.020m 


[□ 


: 


Z 


|0.005m 



Table 4 - Occupations 

Point Name: 
Data file: 
Receiver Type: 
Receiver Serial Number: 
Antenna type: 
Measured To: 

Antenna height Measured 
APC 



Successful Adjustment in 1 iteration(s) 



Network Reference Factor 
Chi Square Test (□ =95%) 
Degrees of Freedom 
GPS Observation Statistics 



1.00 

PASS 

90.00 



3842 

38420790.DAT 

4000SSi 

13842 

Compact L1/L2 
Bottom of antenna 
1.270m 
1.332m 

Table 5 - Adjusted 3D Coordin 
Statistical Summary 



Errors are reported using 1.96 . 



To 

AE000091 
30760790.DAT 
4000SSi 
13076 

Compact L1/L2 
Bottom of antenna 
1.630m 
1.692m 



7°19'02.77547"N 



7°19'20.68657"N 



7°19'20.68709"N 



7°19'20.39312"N 



7°19'19.83339"N 



7°19'18.96313"N 



7°19'17.86287"N 



7°19'17.36933"N 



7°19'17.78599"N 



7°19'20.04550"N 



7°19'21.79386"N 



7°19'22.34130"N 



7°19'22.59900"N 



7°19'23.90880"N 



7°19'25.53640"N 



7°19'27.27521"N 



7°19'28.04179"N 



7°19'29.17996"N 



7°19'30.54718"N 



7°19'31.80050"N 



7°19'33.41452"N 



6°22'00.45644"E 



6°2ril.59110"E 



6 o 21'11.59047"E 



6°21'11.91540"E 



6°21'11.34336"E 



6°21'10.67891"E 



6°21'09.61401"E 



6°21'07.72061"E 



6°21'06.86673"E 



6°21'04.33317"E 



6 o 21'05.17316"E 



6°21'04.59310"E 



6°21'03.57478"E 



6°21'00.67533"E 



6°20'58.39906"E 



6°20'58.24238"E 



6°20'57.19865"E 



6°20'55.42582"E 



6°20'52.51781"E 



6°20'50.60862"E 



6°2049.09372"E 



6°2048.80818"E 



6°2047.82717"E 



6°2048.26109"E 



Northing 

367122.496m 
367666.655m 
367666.671m 
367657.679m 
367640.407m 
367613.581m 
367579.638m 
367564.232m 
367576.927m 
367646.036m 
367699.866m 
367716.615m 
367724.404m 
367764.285m 
367814.010m 
367838.657m 
367867.288m 
367892.524m 
302.460m | 367890.251m 
302.719m| 367899.932m 
302.824m | 367924.792m 
300.627m | 367966.768m 
' 368005.156m 
368054.806m 



Easting 

436929.726m 
435428.264m 
435428.244m 
435438.250m 
435420.771m 
435400.497m 
435367.966m 
435309.940m 
435283.690m 
435205.675m 
435231.224m 
435213.359m 
435182.084m 
435092.965m 
435022.923m 
435018.015m 
434985.875m 
434931.381m 
434842.173m 
434783.560m 
434736.980m 
434728.047m 
434697.792m 
434710.900m 
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AE000024 |7°19'34.79316"N | 6°20'48.74330"E 1 290.426m 368097.230m 434725.519m 
AE000025 |7 0 19'36.66374"N | 6 o 20'46.95549"E 1 273.746m 368154.483m 434670.434m 
AE000026 |7 0 19'37.96679"N | 6°20 , 47.67879"E 1 255.774m 368194.614m 434692.459m 




Figure 1 - Contour of the Site 




Figure 2 - 3D model of the Site 

V.2 Grid Volume Computations 

Volume calculations are made using the average end area method and distances along each roadway alignment. For 
a particular baseline, volume computations may be done for individual roadways, combinations of roadways, or the entire 
design, and may extend through any baseline station limits. Each volume calculation request generates volume calculations 
through the limits specified by the user by accessing the appropriate baseline cross section file to get the end area per station 
(established during the design process) and then determining the volume using the volume 

parameters entered with the design data. A mass ordinate file is created as the volumes are calculated. Any design exceptions 
in the stationing are taken into account in the volume calculations. Corrections for horizontal curvature may be applied to the 
volume computation as desired. 

The determination of earthwork quantities is based upon field crosssections taken in a specified manner before and 
after excavation. Crosssections are vertical profiles taken at right angles to the survey centerline. Every section is an area 
formed by the subgrade, the sideslopes, and the original ground surface. 

Volumes are computed from cross-section measurements by the average end area method. 
\L x (Al + Al)] 

Volume (metric) ( m3) = (16) 

2 

Where L is in meters Al and A2 are in square meters 

Specifications require this calculation. All the plans and bidding for the formula for average end areas is accurate 
only when the end areas are equal. For other cases, the formula generally gives volumes slightly larger than their true values. 
If applied to a pyramid, for example, the error would be the maximum and would be equal to 50 % of the correct volume. 

From fig 2, the minimum upper surface in term of (XYZ) coordinate are (367564.24mN, 434290.1 lmE, 182.542m), 
while the maximum is (368418.331mN, 435438.297mE, 332.052m) at an interval of 12m spacing. The grid size is 100 
rows by 75 columns. The lower surface is defined by Z value of 165m. 

Using Trapezoidal Rule, the total volume of the surface was found to be 128594444.07633m 3 . Using Simpson's 
Rule, the total volume was found to be 128591558.2169m 3 . The positive volume (cut) was found to be 128594039.74377m 3 , 
while the net Volume (cut - Fill) was found to be 128594039.74377m 3 . The total planar area was found to be 
980656.18301705m 2 , while the positive Surface Area (cut) was found to be 996019.8002409m 2 . 
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VI. DISCUSSION AND CONCLUSION 
Our approach present baseline Static and kinematic GPS positioning by using standard Kalman filter process model. 
Table 1 is the baseline processed results. Table 2 is the result of the adjusted 3D coordinate of the site investigated. Figure 1 
is the contour showing the nature of the terrain. Figure 2 is the presented Digital Terrain Model (DTM). The DTM reveal the 
area the will needed to be filled and cut. The quantity of material needed was also computed. The contour, 3D model, the 
area and Volume was derived using SUFFER 8.0 software. 

Performance assessment of the measurement system carried during the data acquisition involved the use of three 
receivers, one control for static model and two rovers for Kinematic model both acquiring data at 5 seconds interval. From 
the results obtained it can be concluded that GPS observation by combining static and kinematic models gives very high 
accuracy for 3-D Survey. 
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